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In “hilltop inflation”, inflation takes place when the inflaton held slowly rolls from close 
to a maximum of its potential (i.e. the “hilltop”) towards its minimum. When the inflaton 
potential is associated with a phase transition, possible topological defects produced during 
this phase transition, such as domain walls, are efficiently diluted during inflation. It is 
typically assumed that they also do not reform after inflation, i.e. that the inflaton field 
stays on its side of the “hill”, finally performing damped oscillations around the minimum 
of the potential. In this paper we study the linear and the non-linear phases of preheating 
after hilltop inflation. We find that the fluctuations of the inflaton field during the tachyonic 
oscillation phase grow strong enough to allow the inflaton field to form regions in position 
space where it crosses “over the top of the hill” towards the “wrong vacuum”. We investigate 
the formation and behaviour of these overshooting regions using lattice simulations: Rather 
than durable domain walls, these regions form oscillon-like structures (i.e. localized bubbles 
that oscillate between the two vacua) which should be included in a careful study of preheating 
in hilltop inflation. 
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1 Introduction 

Inflation has proven to be a very successful paradigm in early universe cosmology. It success¬ 
fully solves the horizon and flatness problems of Big Bang cosmology and at the same time 
provides an explanation for the adiabatic, nearly scale invariant and Gaussian primordial 
curvature perturbations implied by observations of the cosmic microwave background and 
large scale structure [1-5]. 

At the end of inflation, the universe is dominated by the inflaton field’s potential energy 
which is converted into radiation in a process called reheating. Reheating typically happens 
in at least two stages: an initial “preheating” phase of non-perturbative particle production, 
and a later stage of perturbative inflaton decay and thermalization. Together, these processes 
describe the transition of the inflationary, vacuum dominated phase to the radiation dom¬ 
inated universe of standard ACDM cosmology. A detailed understanding of the reheating 
phase is required not only for precise predictions of the primordial spectrum (as the expan¬ 
sion history after inflation affects the matching of inflationary and CMB perturbations), but 
also for calculating the non-thermal production of e.g. a baryon asymmetry, dark matter or 
other relics from inflaton decays. 

Despite its success in cosmology, it is not yet clear how exactly inflation is realized in 
terms of particle physics. One of the many interesting possibilities is hilltop inflation, also 
called “new inflation” (see e.g. [6-11]): the inflaton (j) starts at a local maximum of the scalar 
potential at 0 = 0 and slowly rolls towards a minimum of the scalar potential at 0 ^ 0, 
spontaneously breaking some particle physics symmetry like a family symmetry [12]. The 
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Figure 1. Infiaton potential with the points of maximum tachyonic mass (0m), the end of tachyonic 
preheating (0i) and the potential minimum {v). Inflation happens at the hilltop for |0| < v‘^/ 
as the infiaton field slowly rolls down towards the minimum at (j) = v. 


initial conditions for hilltop inflation (0 starts near a local maximum) can be easily generated 
by a coupling of the infiaton to some matter field x? which allows for a period of preinflation 
[8, 9, 13, 14] and also provides an efficient decay channel for reheating [14]. Other proposals 
for generating the initial conditions include [15-18]. 

Preheating in hilltop inflation has been studied numerically in [19] and analytically in 
[20]. Both papers find that perturbations are amplified by both tachyonic preheating and 
parametric resonance, and that perturbations grow very large at scales corresponding to the 
infiaton mass. Detailed analytical calculations for the linearised perturbations are presented 
in [20]; however, when the perturbations eventually grow too large, the linear approximation 
breaks down and any further calculation must include the non-linear interactions. [19] stud¬ 
ies preheating numerically using lattice simulations, using a slightly different (logarithmic) 
hilltop potential. They study the evolution of particle number densities rik and find that the 
growth of perturbations is terminated by the non-linear interactions when (50^) {4>?- 

In this paper, we study both the linear and the early non-linear stages of preheating in 
hilltop inflation. Beyond confirming the earlier results of [19] and [20], we study the evolution 
of the inhomogeneous infiaton field in position space. We are particularly interested in the 
possibility that the large infiaton fluctuations might push the infiaton field over the local 
maximum of the scalar potential towards the “wrong” vacuum, which could in principle 
result in the formation of domain walls. We find that such “hill crossing” indeed occurs for 
a wide range of parameters. However, it turns out that these regions of “wrong” vacuum do 
not form domain wall networks, but that they instead generate localized oscillating bubbles, 
which might be interpreted as oscillons [21-29]. 

The paper is structured as follows. In section 2, we introduce the hilltop inflation 
model studied in this paper. The linear phase of preheating is discussed in section 3 based on 
[20], with special emphasis on the exponential growth of perturbations which suggests that 
perturbations might grow large enough to push the infiaton field over the hilltop in some 
regions of space. However, the linear approximation breaks down for such large perturbations, 
so to check whether hill crossing can actually happen, we study the non-linear preheating 
dynamics using lattice simulations in section 4. We then summarize our results in section 5. 
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2 Hilltop inflation 

In this paper, we consider the hilltop inflation potential 

y(0) = Fo (i - (2.1) 

with V <C mpi. This potential is plotted in fig. 1. 

Hilltop inflation happens near 0 = 0, with 0 slowly rolling down the hilltop, while the 
universe expands by at least ^ 60 e-folds to solve the horizon and flatness problems. 
We choose 0 > 0, which can always be arranged by a field redefinition 0 ^ —0. Slow-roll 
inflation ends shortly after 77 = /V = — 1, which happens at the inflaton field value 

10^1 ^ and the inflaton quickly rolls towards the minimum at 0 = v.^ 

Primordial perturbations 

The scalar perturbation amplitude matches the observed value Ag = 2.2 x 10“^ for 

l/o^2x 10-^%^ (2.2) 

and the predictions for the spectral index rig and tensor-to-scalar ratio r are 

ns = l-A- 0.95, r ~ 6 X 10"'^ (—) , (2.3) 

\mpij 

in some tension with the most recent Planck bounds rig = 0.965 ± 0.005 (at 68% CL) [1]. 
The predicted value of rig can be larger if there is an additional tachyonic mass term for 0, 
which occurs naturally in supergravity implementations of hilltop inflation due to corrections 
from the Kahler potential, or if the 0^/?;^ in the inflaton potential is replaced by some higher 
power /rP with p > 4. 

Initial conditions from preinflation 

The initial conditions for hilltop inflation (0 ^ 0) can be naturally generated e.g. by a phase 
of preinflation [8, 9, 13, 14]. During preinflation, another scalar field y slowly rolls down 
its potential, while 0 is driven to zero by a large mass term generated from couplings like 
Vint = ^0^X^ [14]. Eventually, when y reaches very small values, the potential for 0 takes the 
form of eq. (2.1), and 0 starts to slow-roll away from zero either towards 0 = 7 ; or 0 = —7;, 
starting with an initial displacement from 0 = 0 due to quantum fluctuations. This second 
slow-roll phase, which can be described as single-field hilltop inflation in 0 with the potential 
(2.1), lasts more than ^ 60 e-folds. As a consequence, any topological defects generated 
during the initial symmetry breaking are inflated away, and at the end of inflation the entire 
observable universe has either 0 > 0 or 0 < 0. 

Preheating 

Preheating in hilltop inflation has been shown to happen in two steps [20]: tachyonic pre¬ 
heating and tachyonic oscillations. 

^In fact, inflation continues for ~ 2.5 e-folds after 0^, but the expansion in potential slow-roll parameters 
breaks down around 77 ~ —1. 
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1. While the inflaton rolls from to (j)\ = (3/7)^^^'?; 0.81?;, the inflaton mass is tachy- 

onic: V'\(j)) < 0. During this phase, long-wavelength perturbations with A: < i7o grow 
exponentially in a process called “tachyonic preheating”. 

2. For V > 10“^mpi, (j) does not quickly relax at the minimum. Instead, it oscillates 

between 0^0 and 0 ^ These “tachyonic oscillations” cause a strong growth in 

small-wavelength perturbations around a scale k ^ /Cpeak — 7(mpi/?;)^/^i7o, which grow 
until the perturbations become non-linear (for v/m^i < 10“^) or until Hubble damping 
reduces the oscillation amplitude so much that the inflaton oscillations remain in the 
approximately quadratic part around the minimum. 

3 Linear preheating 

In this section, we discuss the phase of linear preheating after inflation. During this phase, 
the inflaton field is still almost homogeneous with small fluctuations 5(/) which can be treated 
as small perturbations (keeping only terms to leading order in Sc/)). The initial phase of 
preheating can be studied by solving the time evolution of these fluctuations in Fourier 
space. 

We start with a brief review of the theory of linear preheating and its application to the 
hilltop inflation potentials in eq. (2.1). For these potentials, preheating starts with tachyonic 
preheating, followed by a period of tachyonic oscillations. We summarize the analytic results 
from [20] for these two phases and show numerical results for the spectra produced by linear 
preheating. We then discuss how these power spectra suggest that for 10“^ < ^ 10“^, 

the perturbations could push the inflaton field above the hilltop towards (/){x) < 0, a process 
which could potentially produce domain walls. 

3.1 Introduction to linear preheating 

3.1.1 Equations of motion 

The equation of motion of the inflaton field (/)(t, x) in a FRW background (neglecting back- 
reaction from the metric perturbations^) is given by 

^ - VV + V\(t)) = 0. (3.1) 

At the end of inflation, the inflaton field is almost homogeneous, so we can expand it in terms 
of a homogeneous mean field and a perturbation field 5(j){t^x)\ 

(j){t^x) = (j)(t) + 5(j)(t^x). (3.2) 

As long as the perturbations remain small, we can linearise eq. (3.1) by dropping all terms 
of 0(502): 


0 + 3770 +1/'(0) = 0, (3.3a) 

^To account for scalar metric perturbations, one can use the Mukhanov-Sasaki variable instead of the 
inflaton held value [30, 31]. For hilltop inhation, the effect on the tachyonic preheating spectrum is known 
to be negligible [20]. We have also checked for a few parameter points that including metric perturbations in 
this way did not make a noticeable difference for the subsequent oscillation phase, so analogously to [20], we 
perform our calculations using the scalar held eq. (3.1), neglecting metric perturbations. 
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(3.3b) 


54> + 3H5<j) + (v"{(p) - 5(p = 0 . 

The partial differential equation (3.3b) for 5(j){t^x) can be transformed into ordinary differ¬ 
ential equations by transforming to Fourier space: 

4>k + 3H<Pk + {v"{^) + ^) = 0, (3.4) 

where (j)]^ is the Fourier mode with comoving momentum k. 

3.1.2 Initial conditions from Bunch-Davies vacuum 

The initial conditions are given by matching to the Bunch-Davies vacuum at early times: 

lim (pk{t) — / (i (^3 5 ) 

For a numerical implementation, where matching at t ^ —oo is not practical, “early times” 
means that the Hubble constant at the time of matching must still be approximately constant, 
and that the inflaton mass V"{(j)) should be negligible compared to k‘^ jo?. 

3.1.3 Inflaton fluctuation amplitude and power spectrum 

Given eqs. (3.4) and (3.5), one can calculate the evolution of the scalar field perturbations 
for each mode k. The result can be translated into position space perturbations by an inverse 
Fourier transformation. For the inflaton field’s variance, this gives 

((5^^(t,f)) = J {dlnk)V^{t,k), (3.6) 

where the limits of the integral are the IR and UV cutoffs,^ and for d = 3 spatial dimensions 
the spectrum can be calculated as 

'P4>i^,k) = . (3.7) 

The results from linear preheating are reliable only until higher orders in 50 become 

important. At that point, the fully non-linear eq. (3.1) must be solved. This is usually done 

by solving the equation numerically on a discrete spacetime lattice, as we will discuss in 

section 4. For the hilltop potential (2.1), we can estimate that the linear approximation is 

valid approximately until (50^) > 0^, or equivalently maxP^ > cjp. 

k 

3.2 Tachyonic preheating 

When R'X0) < 0, all modes with k/a < y/—'F"(0) experience exponential growth. This 
period of tachyonic growth lasts the longest for modes with A; < ido, which are already 
tachyonic at the end of inflation. Higher k modes become tachyonic later, and therefore 

®The infrared cutoff can be naturally introduced by considering a box of finite size, e.g. a single Hubble 
patch or the observable universe. The ultraviolet cutoff must be put in by hand to regularize the usual 
UV divergence arising from small-scale quantum fluctuations. This cutoff should be chosen large enough to 
encompass the largest (pk amplihed by preheating dynamics, but small enough for the non-amplihed quantum 
fluctuations to be subdominant. 
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Figure 2. Oscillations of the inflaton background field around the minimum as a function of Hot, 
plotted for v/rripi = 10“^ and vjm^x = 10“^. The dashed horizontal line denotes the inflaton value 
O.8I1; below which V"{(j)) < 0. We see that for small v/mp\ 10“^, the oscillations continue to 
bring 0 deeply into the tachyonic region for several oscillations, which explains why the perturbations 
grow non-linear within a few oscillations. For large v/rripi = 10“^, however, Hubble damping quickly 
dampens the oscillation amplitude so that after a few oscillations, the inflaton field never re-enters 
the tachyonic region. Preheating is much less efficient in this case: the perturbations remain small 
and are eventually redshifted away as the inflaton field oscillates around the quadratic part of its 
potential. 


achieve a smaller final amplitude because they grow for a shorter time. Modes above k^/a = 
max y/—!/"((/)) 4:{mp\/v)Ho are not amplified at all. 

After tachyonic preheating^, the power spectrum around the Hubble scale is 

V^iaHo) - 10-^^ mix, (3.8) 

independent of v. For aHo < A: < /cm, the power spectrum falls off as oc {k/a)~‘^, and for 
k ^ fcm, the spectrum is just the vacuum spectrum /c^/(47r^a^). 

The variance of the inflaton field is dominated by the large contribution from the Hubble 

scale: 


(Scp^ix)) = j {d\ogk)V^{k) -- V(i,{aHo) -- 10 (3.9) 

For very small v/rripx % 10“^, this implies {dcjp^x)) > so the perturbations become 
non-linear during tachyonic preheating - the inflaton field becomes very inhomogeneous and 
the description in terms of a background field plus perturbations breaks down. For larger 
^ 10“^, eq. (3.9) implies (50^) <C and the linear approximation should remain 

viable. 

3.3 Tachyonic oscillations 

During the subsequent phase of tachyonic oscillations, the background field oscillates around 
the minimum. For vjvripx ^ 10“^, these oscillations take the inflaton back into the tachyonic 

^The value of the power spectrum in eq. (3.8) is evaluated at 0 = u at which time the tachyonic peak takes 
on its maximum. 
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region with 0 < 01. In this case, one can estimate the inflaton field value after j complete 
oscillations as 


(3.10) 

V \ mpiy 

These oscillations trigger strong growth of perturbations at a characteristic peak scale 

c, (3.11) 

For small v/m^\ <C 10“^, the perturbations grow by many orders of magnitude and become 
non-linear within a few oscillations around the minimum. For large v/m^\ > 10“^, however, 
the oscillations are quickly damped due to Hubble damping (see fig. 2). After a few oscilla¬ 
tions, the amplitude decreases so much that the inflaton field never re-enters the tachyonic 
region, which makes the preheating much less efficient. In this case, the perturbations re¬ 
main so small that we expect the linear approximation to remain valid throughout the entire 
preheating phase. ^ 

Note that the perturbations with k <C A^peak? which are amplified most during tachyonic 
preheating, do not continue to grow during the oscillation phase. Even though they are 
amplified every time that the inflaton field rolls down the tachyonic part of the potential, 
this amplification is exactly cancelled by an exponential damping sustained every time the 
field rolls back up. The infrared part of the spectrum therefore oscillates: it is very large 
near the minimum at (p ^ and it is negligible at pj (i.e. near the hilltop). 


3.4 Numerical results for the power spectra 

Integrating the mode equations (3.4) numerically, we can calculate the power spectrum ^^^(fc) 
of the inflaton perturbations. Fig. 3 shows for various values of v evaluated at the 

time when the inflaton crosses the minimum (p = v and 0 > 0), while fig. 4 shows the power 
spectra closest to the hilltop at the end of each full oscillation, where 0 = 0. 

All of the plots show the strong growth of power for modes with k ^ /Cpeak with each 
successive oscillation. For v/m^i < 10“^, it only takes a few oscillations for the perturbations 
to grow non-linear (i.e., ^^^(/Cpeak) ^ 0^)- Only for large v/mp\ — 10“^, the perturbations 
remain small due to Hubble damping, as explained above. 

For the infrared modes, fig. 3 clearly shows the tachyonic preheating spectrum with 
V(f){aHo) ^ and V(i){k) oc k~^. The infrared modes do not continue to grow during 

the tachyonic oscillation phase: for v/m^\ < 10“^, the spectra after the j-th oscillation lie 
on top of each other. For v/m^\ — 10“^, we even see that the infrared modes slowly decay 
due to Hubble damping and redshifting. Near the hilltop, the tachyonic preheating spectrum 
almost vanishes, as one can see in fig. 4. The infrared part of the spectrum actually oscillates 
during the tachyonic oscillation phase, growing to maximum values at 0 ^ 'z; and shrinking 
to minimal values at 0 0j, whereas 'P(/)(A:peak)/0^ continuously grows throughout each 

oscillation^ until the evolution becomes non-linear. 

® There is some non-adiabatic growth of the inflaton perturbations during the subsequent non-tachyonic 
oscillations around the minimum, but red-shifting and Hubble damping make the perturbations decay before 
they can grow non-linear. 

^To be precise, ^^(/cpeak)/^^ temporarily drops when cj) > v, but it grows both on its way up and down 
the hilltop for (j) < v. 
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Figure 3. Linear power spectra evaluated at (p = v and 0 > 0 for the first few oscillations (first to 
last oscillation: blue to red). Tachyonic preheating enhances modes with k < /Cmax — O.dSSmmin, with 
the strongest amplification in the infrared {k < Hq), whereas the oscillations enhance modes with 
much larger momenta around the peak scale /cpeak — which is denoted by a vertical 

dashed line. For v/rripi = 10“^, we do not show every oscillation, because linear preheating lasts for 
very many oscillations; the dark red line in this case corresponds to over 1000 oscillations. 


3.5 End of linear preheating and hill crossing 

Linear preheating predicts that for 10“^ < v/'ni.pi ^ 10“^, there is a rapid growth of pertur¬ 
bations at the scale fcpeak? producing a large fluctuation amplitude 

{S4>‘^{x)) = J{dlnk)V4,{k) -- P,^(/cpeak) ^ (3.12) 

This could indicate either that the fluctuations become large enough to push the inflaton 
field over the hilltop towards (/){x) < 0, or that linear perturbation theory breaks down and 
non-linear interactions stop the growth of perturbations before this happens. To find out 
which of these options is realized, we need to go beyond the linear approximation. 

As we will show below, it turns out that the fluctuations can really push (/)(x) to negative 
values at some points in space, and the potential subsequently drives (/){x) —v in these 
regions. The inflaton field is then no longer well described by a homogeneous field 0 plus small 
perturbations 5(/){x). To study the formation and evolution of the regions with (/){x) < 0, 
we need to solve the field equations of motion non-perturbatively using numerical lattice 
simulations. 
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Figure 4. Linear power spectra normalized to the mean infiaton value evaluated at the end of each 
oscillation when the infiaton is closest to the hilltop (first to last oscillation: blue to red). We see that 
the infrared peak from tachyonic preheating is negligible near the hilltop, whereas the perturbations 
around the peak scale /Cpeak (vertical dashed line) can grow amplitudes larger than cpj , indicating that 
on these scales, the perturbations can push the infiaton towards (j){x) < 0. For v/m^\ = 10“^, we 
do not show every oscillation, because linear preheating lasts for very many oscillations; the dark red 
line in this case corresponds to over 1000 oscillations. 


4 Non-linear dynamics 

In this section, we present and discuss the results of numerical lattice simulations of the 
non-linear stage of preheating for v/m^\ — 10 “^ and v/m^\ — 10 “^. 

In the previous section, we have seen that for 10“^ < v/m^\ < 10“^, the infiaton 
fluctuations get amplified during the initial descent towards (j) = v dX scales fctac ^ aiLo- 
Then the system oscillates between cj) v and (j) > v during which time the spectrum 
gets amplified at scales /Cpeak ^ either until non-linearities become dominant 

(for v/m^\ < 10“^) or until the oscillations become damped and redshifted by the Hubble 
expansion (for v/m^\ > 10“^). 

In order to study the non-linear phase of preheating, we numerically solve the non-linear 
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equations of motion 


0(t, x) + x) — \v^(/)(t, x) + ^ = 0 , (4.1) 

0(p 

where V are gradients with respect to the comoving coordinates x. We use a modified version 
of the program LATTICEEASY [32] to solve these equations on a discrete spacetime lattice. 

The hierarchy between /ctac sind /Cpeak implies that the lattice must have much more 
than points per dimension to include both scales in the simulation. For small 

v/m^\ 10 this becomes impractical due to our limited computing power. 

However, from the linear analysis we expect the fluctuations at fctac fo remain relatively 
small, as only the fluctuations at the peak scale /Cpeak grow during the oscillations in the linear 
phase. At the end of the linear phase, most of the perturbations will be around k ^ fcpeak- 
We therefore ignore the large wavelength modes at /ctac and only include scales around /cpeak 
in order to resolve the dominant small distance fluctuations. 

For v/m^i = 10“^, we could even keep most of the infrared modes, because the hierarchy 
between /ctac and /Cpeak is smaller. However, the linear power spectrum at ktac is very small 
for v/m^i — 10“^, so even in this case we prefer to cut off some infrared power in favour of a 
larger UV cutoff, which is more important to resolve the small-distance fluctuations produced 
by the hill crossing. 

To resolve a reasonable range around /Cpeak? we focus on 2 spatial dimensions, although 
we also did some lower-resolution runs in 3 dimensions (see section 4.3) and saw the same 
qualitative behaviour. 


4.1 Initialization of the lattice 

The initial time of the lattice simulation must be chosen such that the system is still well 
described by the linear equations (3.4). We choose to start after r] = —1 but shortly before 
6: = ^^pi(EYE)^ = 1, with initial conditions and parameters shown in table 1. 

The initial field fluctuations (0^)iin for = 10“^mpi and v = 10“^mpi are extracted 
from the linear equations (3.4). As shown in fig. 5, the power spectra around fcpeak 
time are mainly dominated by the vacuum contribution. The tachyonic contribution becomes 
important for scales k < O.l/Cpeak? corresponding to the infrared tail of our initial spectra. 

As usual, the norm of the fluctuations is initialized as a stochastic variable obeying the 
probability distribution [33, 34] 


PD{\4>k\) = 


2|0fc| 

{<Plhin 


exp 



(4.3) 


The initial Fourier components of 0 and their derivatives are given by 


cf>k = t + , 

4 = — , - H^k ■ (4.4) 
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v=10 ^mpi 


v=10 ^mpi 




Figure 5. Initial spectra of our lattice simulations for v — 10“^mpi (left) and v — 10“^mpi 

(right), with other parameters given in table 1. The vertical dashed line corresponds to k = /cpeak 
and the range of scales goes from k[r = /Cuv/4096 to k^v = 100/cpeak- For k > O.l/cpeak, the spectra 
are dominated by vacuum fluctuations T^vac oc whereas for k < 0 .1/cpeak one can see the effect of 
tachyonic preheating. The scale /ctac is outside of the range for both v = 10“^mpi and v = 10“^mpi, 
which we expect not to affect our result significantly as perturbations on that scale are subdominant 
compared to fluctuations at /Cpeak (see section 3). 


vlm^i 


(0)iY 

^UV 

N 

LHo 

10“2 

0.058 

8.9 X 10“^° 

100/Cpeak 

4096 

0.83 

10-5 

0.0098 

2.7 X 10-^1 

100/Cpeak 

4096 

4.7 X 10-5 


Table 1. Initial field and derivative mean values and other parameters of the 2D lattice simulations. 

where cr_, 6^ and 9- are random real numbers uniformly distributed between 0 and 1.^^ 
In setting the derivative we assumed that (^k oc Although this is true only 

for the vacuum fluctuations and not for the infrared part k < O.l/Cpeak of fhe the spectra in 
fig. 5, our results are not affected by this assumption. Indeed, the lattice simulations’ results 
during the linear phase match those of the linear analysis. 

The fluctuations (4.4) are then Fourier transformed to position space and a discretized 
version of the equations of motion is solved on a lattice with spacing Sx = with 

fcuv = lOO/Cpeak as the ultra-violct cutoff. The size of the lattice is given by L = jk^^ 
where N is the number of points per dimension. This defines the infrared cutoff k^^ — k^-^/N. 
In section 4.2 we present results of simulations in T) = 2 spatial dimensions with N = 4096 
and in section 4.3 we briefly discuss simulations in T) = 3 spatial dimensions with N — 256. 

In order to run the simulations, we must choose suitable rescalings. We define our 
program units to be 


0pr — 0/^0, Xpr — ^peak^ 5 dtpr — ^peak^^ 5 (^•^) 

This corresponds to choosing LATTICEEASY rescalings A — B — /Cpeak? r — 12 and 

5 = 0. 

^°Note that this is different from the original LATTICEEASY initialization, where the fluctuations are 
initialized as standing waves, corresponding to a- — = 1. 
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4.2 Results of lattice simulations in 2 dimensions 

Based on the linear analysis of section 3.4, we expect that for v = 10“^mpi the dynamics 
remain linear throughout the oscillations phase. Our lattice simulations confirm this result, 
with the variance growing to a maximal value of {dcjP) ~ . 

For V < 10“^mpi, on the other hand, preheating is expected to become non-linear already 
during the tachyonic preheating stage [20], which should prevent a tachyonic oscillations 
phase, so that no spectral peak at k^eak is expected. 

For these reasons, we focus on the lattice simulation results for v/mp\ = 10“^ and 
v/m^i = 10“^. In 2 spatial dimensions we can include N = 4096 points per dimension, which 
is sufficient to resolve the relevant scales around /Cpeak foi* both of our choices of v. 

4.2.1 Position space slices 

Figs. 6-7 show the field values on the 2D lattice at 9 different points in time, for v = 10“^mpi 
and V = 10“^mpi, starting around the time when the initial hill crossing occurs.The red- 
coloured regions have field values around the positive minimum (j) ^ the yellow and light 
blue regions are close to the hilltop 0^0, and the dark blue regions are around the negative 
minimum (j) ^ —v. Animated movies showing many more time steps for both simulations are 
available at [35]. 

We clearly see that overshooting to 0 < 0 indeed occurs, creating small negative regions 
with 0 < 0 which are separated by a distance of roughly Apeak = ‘^^/kpeak- The size d- of 
these negative regions is significantly smaller than the typical distance Apeak? especially for 
V = 10“^mpi. This is what one would expect from overshooting due to fluctuations dominated 
by a specific wavelength Apeak- fhe regions which overshoot are those close to the maxima of 
this wave, which explains both their small size d- < Apeak ^^nd the typical distance Apeak- 

During the early stage around the initial overshooting, the dynamics is very different 
for V = 10“^mpi and v = 10“^mpi. For v = 10“^mpi, the initial overshooting occurs while 
(0) ^ V. The initial negative regions are small, and they very quickly start oscillating towards 
(j) —V. For V = 10“^mpi, the initial overshooting occurs close to the hilltop with (0) <C 

with larger initial size of the negative regions. Although the first three steps in fig. 7 look 
similar to the early stages of the formation of a domain wall network, the subsequent evolution 
is qualitatively different (best seen in the animated movie [35]) and instead leads to many 
isolated bubbles with <p ^ —v completely enveloped in a single connected region with <p ^ v. 

In both cases, we eventually end up with small overshooting bubbles (blue), roughly 
separated by the distance Apeak? which are connected by filaments in which the field is close 
to the hilltop (yellow). These filaments then vanish, and only the blue overshooting bubbles 
remain. 

For V = 10“^mpi, we can clearly see that the overshooting bubbles oscillate between 
0 < 0 and 0 > 0. These oscillations are localized and mostly in phase: we can clearly see 
that the blue regions vanish in the 5th, 7th and 9th lattice slice, and most of them reappear 
in the 6th and 8th slice at approximately their original positions. Over time, the oscillations 
get damped and some of the blue regions evaporate. The oscillations of the different bubbles 
also get out of phase over time, which is why the 9th slice has more blue bubbles than the 
5th: in the 5th slice, all bubbles have just oscillated towards 0 > 0 at the same time, whereas 

^^Note that, to minimize image sizes, the slices shown in figs. 6-7 are from a simulation with N = 1024, 
compared to N = 4096 of figs. 8-9. Nevertheless, the initial fluctuations are the same in both simulations and 
the power lost in the infrared is subdominant to that at the peak scale /cpeak, as one can see from fig. 9. 
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in the 9th slice, some bubbles are lagging behind (or ahead) and are still (or already) at 

For V — 10“^mpi, the blue bubbles in fig. 7 do not clearly exhibit any coherent oscil¬ 
lations. The largest bubbles indeed tend to not oscillate but quickly fragment into smaller 
structures. However, some of the smaller blue bubbles oscillate between 0 < 0 and 0 > 0, 
but these oscillations are strongly out of phase so that fig. 7 mostly shows the overall decay 
of the oscillations.^^ 

We might interpret these oscillating bubbles as “oscillons”. It has been shown in the 
literature that such oscillons can form from collapsing bubbles in double-well potentials, and 
that they can be very long-lived, depending on their properties (like radius and initial field 
displacement) [21]. While we observe that some of the oscillon-like bubbles with 0 < 0 vanish 
within a few oscillations, our simulation does not extend to sufficiently late times to check 
whether the others might remain stable over a longer timescale, or whether the bubbles 
which disappeared continue as long-lived oscillons which are not “hill crossing”. We note 
that recently, for a double-well potential similar to ours, oscillons (not “hill crossing” ones) 
have been observed which remain long-lived on cosmological timescales [28]. 

In summary, we find that hill crossing occurs and produces regions with 0 < 0. However, 
these regions do not form domains settling dX (j) = —v separated from the other vacuum by 
domain walls, but instead they behave like localized oscillations between (j) ^ —v and (j) ^ v. 
While some of these oscillon-like bubbles evaporate over time, our simulation does not clearly 
indicate whether or not some of them may be longer-lived. 

4.2.2 Mean, variance and fraction of negative field regions 

The left plots in fig. 8 show the mean (0) and perturbation amplitude of the infla- 

ton field. We see the first few coherent oscillations of the background inflaton field during 
which the perturbations remain relatively small. When the perturbation amplitude becomes 
comparable to the mean, the coherent background oscillations break down and hill crossing 
occurs. This happens during the 7th oscillation for v — 10“^mpi and during the 3rd oscil¬ 
lation for V — 10“^mpi, in agreement with our expectations from the linear analysis. For 
V — 10“^mpi, some coherent oscillations in (0) and {dcjP) remain after the hill crossing, as 
the negative bubbles oscillate in phase between 0 < 0 and 0 > 0. For v — 10“^mpi, the 
averages (0) and (50^) do not oscillate because the individual bubbles’ oscillations are out 
of phase and therefore roughly average out. 

The hill crossing is most obvious in the right plots of fig. 8, which show the ratio of 
points with 0 < 0 as a percentage of total points on the lattice. For v = 10“^mpi, this plot 
also clearly shows both the in-phase oscillations of these points between 0 < 0 and 0 > 0, 
and the overall decay of the negative regions over time. One can also see that the oscillations 
gradually get out of phase: this is why the minima in this figure become larger during later 
oscillations. 

For V = 10“^mpi, the right plot of fig. 8 shows no clear oscillations because the individual 
bubbles oscillate out of phase. However, one can still see both the initial overshooting and 
the subsequent decay of the overshooting bubbles. 

We also see that just after hill crossing, the ratio of points with 0 < 0 is larger for 
smaller v. We have observed this general tendency in several lattice runs, though the exact 
ratio varies with the random phases chosen when initializing the lattice, see section 4.1. 

^^The oscillations are clearly visible in the animations containing more time slices, which are available online 
[35]. 
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4.2.3 Power spectrum 

Fig. 9 depicts the power spectra for different points in time, starting from around the time 
of hill crossing. We have checked that until shortly before hill crossing, the power spectra in 
the lattice simulation match the results from the linear analysis. 

Around the time of hill crossing, we can see that the non-linearities wash out the oscilla¬ 
tory features in the spectrum for k > /Cpeak? sind that they transfer power from /Cpeak fo larger 
k. For example, hill crossing transfers some power by producing small overshooting regions 
of size d- < Apeak from perturbations of a longer wavelength Apeak- As explained above, 
this is expected because overshooting should happen near the local maxima of the dominant 
perturbation modes, and such regions around the local maxima are necessarily smaller than 
the wavelength. 

For V = 10“^mpi, the spectrum remains peaked at /Cpeak? whereas for v = 10“^mpi, the 
spectrum develops much more power at small scales k > /Cpeak- This fits well to the lattice 
slices in figs. 6-7, where for v = 10“^mpi the bubbles roughly retain the size they had during 
the initial overshooting, whereas for v — the bubbles tend to fragment and smaller 

features develop. 

Eventually, the spectrum flattens over the /c-range covered by the lattice (especially for 

V = 10“^mpi). This limits the maximum time for which our lattice calculations remain valid. 
Extending the lattice simulation beyond that time would require more lattice points and 
therefore much more computing power. 

4.3 Results of lattice simulations in 3 dimensions 

In order to confirm the validity of our results in 3D, we repeated the simulations for v = 
10“^mpi and v = 10“^mpi with N — 256^ points. Because of the lower resolution, we 
decreased the ultraviolet cutoff to k^y = 50/Cpeak- The infrared cutoff is therefore kir ^ 
0.2/cpeak- The initial conditions for the field mean and fluctuations are the same as in the 2D 
simulations, see table 1 and fig. 5. 

The 3D simulations show the same qualitative behaviour as the 2D simulations. Hill 
crossing happens shortly after the end of the linear phase, and we reproduce fig. 8 throughout 
the coherent oscillations and the initial overshooting phases. For v = 10“^mpi, the 2D 
results are reproduced until the end of the simulation at t = ^^^peak' other hand, 

for V = 10 “^mpi, the simulation is only reliable until t ^ ^^^peak^ because afterwards the 
spectrum is dominated by the smallest scale fluctuations in the lattice. 

Time slices of the 3D simulations are shown in fig. 10.^^ For both v/mp\ = 10“^ and 
10“^, the initial hill crossing phase is characterized by filament-shaped regions where 0 < 0. 
Eventually the filaments disappear and are replaced by localized oscillations between (j) ^ —v 
and (j) with typical separation Apeak = 27r/fcpeak5 as we saw in the 2D simulations. 

Furthermore, for v — 10“^mpi, we can see a tendency of most overshooting regions to 
develop towards a spherical shape. For v = 10“^mpi, the overshooting regions are too small 
for their shape to be accurately resolved in the lattice in 3D. 

Recapitulating, the 3D simulations confirm the results of the 2D simulations. The linear 
and initial hill crossing phases are accurately reproduced for both v/mp\ — 10“^ and 10“^. 
For V — 10“^mpi the 3D simulation remains well-resolved until the final time of the 2D 
simulation and we see the overshooting regions form and tend towards a spherical shape. For 

V — 10“^mpi, the resolution becomes insufficient shortly after the initial hill crossing. 

animation of the 3D simulation for v — 10“^mpi showing more time slices is available online [35]. 
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5 Summary and conclusions 


In this paper, we have studied preheating after “hilltop inflation”, where inflation takes 
place when the inflaton field rolls slowly from close to a maximum of its potential (i.e. 
the “hilltop”) towards its minimum. When the inflaton potential is associated with a phase 
transition, possible topological defects produced during this phase transition, such as domain 
walls, are efficiently diluted during inflation. It is typically assumed that they also do not 
reform after inflation, i.e. that the inflaton field stays on its side of the “hill” (where without 
loss of generality 0 > 0), finally performing damped oscillations around the minimum of the 
potential. 

To investigate the preheating stage, we first integrated the mode equations during the 
linear phase of preheating, during which inhomogeneities are still small enough to be described 
as small perturbations on a homogeneous background. When the inhomogeneities get larger, 
i.e. during the subsequent non-linear stage, we solved the inflaton field’s classical equations 
of motion with lattice simulations in 2D and 3D. We were particularly interested in the 
possibility that large inflaton fluctuations during preheating might push the inflaton field 
over the local maximum of the scalar potential towards the “wrong” vacuum, which could in 
principle result in the formation of domain walls. 

We found that for 10“^ < < 10“^ the fluctuations of the inflaton field during 

the tachyonic oscillation phase indeed grow strong enough to allow the inflaton field to form 
regions in position space where it crosses “over the top of the hill” towards the “wrong 
vacuum” where 0 < 0. The negative regions with 0 < 0 are separated by a distance of 
roughly Apeak = 27r/A;peak5 one would expect from “hill crossing” due to fluctuations 
dominated by a specific wavelength Apeak- Rather than forming durable domain walls, these 
regions develop into localized bubbles that oscillate between (j) ^ —v and (j) ^ which might 
be interpreted as oscillons [21-29]. These oscillon-like bubbles should be included in a careful 
study of preheating in hilltop inflation. 

As next steps beyond the work presented here, it will be interesting to study the evolu¬ 
tion of these oscillon-like structures in more detail, to study their fate (and lifetime) after the 
hill crossing phase. Furthermore, it will be interesting to study the cosmological consequences 
of the hill crossing phase, especially when a second “matter field” is included which couples 
to the inflaton and receives its mass dynamically from the vacuum expectation value of the 
inflaton field (as e.g. in [14]). This may lead to explosive matter particle production at the 
hill crossing regions, but may also trigger a faster decay of the oszillating bubbles. Finally, 
in order to derive precise predictions of explicit hilltop inflation models it will be crucial to 
reliably calculate the expansion history of the universe, which depends on the evolution and 
lifetime of the hill crossing regions. 
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'L’ = 10 ^mpi, N = 1024, D = 2^ = 100/cpeak- 




Figure 6. (j){x) from the lattice simulation at times 44.25, 44.75, 45.25, 46, 47.1, 49.23, 51.1, 53.5, 
55.2 in units of The corresponding ratios of points with 0 < 0 can be read off fig. 8, where the 

slices correspond to the blue dots in the upper right plot. The legend is shown in units of v. The 
simulation clearly shows the initial overshooting, formation and dissipation of a filament structure, 
and eventually localized bubbles oscillating in phase between 0 ^ —v and 0 ^ u, as discussed in 
section 4.2.1. A higher time resolution movie of the evolution can be downloaded at [35]. 
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Figure 7. 0(x) from the lattice simulation at times 25.55, 26, 27.1, 28.1, 29.55, 31.56, 33, 35, 37 
in units of The corresponding ratios of points with 0 < 0 can be read off fig. 8, where the 

slices correspond to the blue dots in the lower right plot. The legend is shown in units of v. The 
simulation clearly shows the initial overshooting, formation and dissipation of a filament structure, 
and eventually localized bubbles with 0 < 0 which decay over time, as discussed in section 4.2.1. 
Many bubbles also oscillate between (j) ^ —v and (/) ^ v (especially the smaller ones), though their 
oscillation is out of phase and thus only visible with a higher time resolution. A higher time resolution 
movie of the evolution can be downloaded at [35] . 
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'L’ = 10 ^mpi, N = 4096, D = 2^ k^y = lOOfcpeak- 




Figure 8. Left: The variance ^/{S^) and the mean (0) in units of as a function of /Cpeak^ foi* 
V = 10“^mpi (above) and v = 10“^mpi (below). We can clearly see how the mean performs six or two 
oscillations, respectively, before the variance grows to values {Sc/)) ^ (0) and non-linearities become 
important. This matches the expectations from the linear analysis, see figs. 3 and 4. 

Right: The fraction of lattice points where 0(x) < 0 as a function of /cpeak^ for 'T' = 10“^mpi (above) 
and V = 10“^mpi (below). The blue dots correspond to the slices in figs. 6 and 7. 

We clearly see that around the onset of the non-linear stage, hill crossing happens for both values of 
V. For V = 10“^mpi, the crossing regions oscillate (initially in phase) between (f ^ —v and (f ^ v, 
which makes the ratio of points with 0 < 0 oscillate. Over time, the amplitude of this oscillation is 
damped, as some of these oscillating bubbles decay. They also get out of phase, which explains why 
the minima of the ratio of negative points get larger during the later oscillations. For v = 10“^mpi, the 
oscillations are less pronounced: some of the smaller bubbles oscillate, larger bubbles tend to fragment 
into smaller substructures. In this case, the bubbles’ oscillations are also out of phase already initially, 
which is why the oscillations are not clearly visible in the plots (which average over all of the many 
bubbles in our lattice). 
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V = 10“^mpi, N = 4096, 4? = 2, k^y = 100/cpeak- 



Figure 9. The spectrum 7^^, defined by ((^0^) = /(din/c)7^0. 

For V = 10“^mpi (above), the spectra are shown at times 40, 44, 47 and 60 (blue to red) in units of 
^peak* spectra are progressively shifted towards the ultraviolet, although the tail is still dominated 
by the vacuum fluctuations at the end of the simulation. We conclude that the relevant scales are 
well resolved throughout the entire simulation. 

For V = 10“^mpi, the spectra are shown at times 25, 26, 28, 30, 32 and 40 (blue to red) in units 
of The ultraviolet tail progressively grows and by the end of the simulation the peak of the 

spectrum is only one order of magnitude larger than the ultraviolet tail, which limits our ability to 
study the late-time behaviour of the overshooting bubbles. However, during most of the evolution the 
neglected UV power is much smaller than the power at the well-resolved scales around /cpeak- 
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V = 10 ^mpi, N — 256, D — 3, k^v = 50fcpeak- 


V = 10 ^mpi, N = 256, D = 3, kuv = 50fc] 


peak • 


Figure 10 . 0(x) at fixed times from the 3D lattice simulations for v = 10“^mpi (above) and v = 
10“^mpi (below). The contours correspond to points where 0 = 0 (yellow) and 0 = —0.8u (blue). 

For V = 10“^mpi, the left box corresponds to t = 45.and 5.8% of points with 0 < 0, and 
the right box to t = 53.and 3.2% of points with 0 < 0. The initial filament-shaped negative 
regions evolve into bubbles with a tendency towards a spherical shape. As in the 2D simulations, 
these bubbles perform localized oscillations between 0 ^ —v and 0 ^ u. The oscillations are visible 
in the animation containing more time slices, available online [35]. 

For V = 10“^mpi, the left box corresponds to t = and 10.7% of points with 0 < 0, and the 

right box to t = 30/^“^!^^^ and 4.4% of points with 0 < 0. The initial filament-shaped negative regions 
are larger than for v = 10“^mpi. They fragment to smaller structures in a manner similar to our 2D 
simulations; the right box matches the corresponding 2D result, which is roughly the middle slice in 
fig. 7. However, due to the limited resolution in 3D we cannot resolve the subsequent evolution well 
enough to ascertain that they form oscillating bubbles as they did in 2D. 
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